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ABSTRACT 

High-resolution data of the peculiar magnetic massive star HD 148937 were obtained with 
Chandra-HETGS, and are presented here in combination with a re-analysis of the older XMM- 
RGS data. The lines of the high-Z elements (Mg, Si, S) were found to be unshifted and relatively 
narrow (FWHM of about 800kms~^), i.e. narrower than the O line recorded by RGS, which 
possibly indicates that the hot plasma is multi-thermal and has several origins. These data further 
indicate a main plasma temperature of about 0.6 keV and a formation of the X-ray emission at 
about one stellar radius above the photosphere. From the spectral fits and the H-to-He line ratios, 
the presence of very hot plasma is however confirmed, though with a smaller relative strength than 
for the prototype magnetic oblique rotator 6^ OriC. Both stars thus share many similarities, but 
HD 148937 appears less extreme than 0^ Ori C despite having also a large magnetic confinement 
parameter. 

Subject headings: X-rays: stars - stars: massive - stars: early-type - stars: individual (HD 148937) 



1. Introduction 

The peculiar class of Of?p stars has received 
much attention in recent years. These objects were 
first distinguished from normal Of supergiants 
some 40 years ago (IWalbornlll972h because of the 
unusual presence of CiiiA4650A in emission with 
a strength comparable to the neighbouring Niii 
lines. In the last decade, additional peculiar prop- 
erties were discovered such as recurrent spectral 
variations (in the Balmer, Hei, Cm, Sim lines), 
large X-ray over luminosities, and the presence of 
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strong magnetic fields, leading to a refinement in 
the definition of the class (see Naze et aLll2008lj : 
Walborn et al. 2010l . and references therein). In 
total, five Of?p stars are now identified in the 
Galaxy. 

HP 148937 was oii e of the first two Of?p iden- 
tified (IWalbornlll97l . but it was only studied in 
detail recently. Contrary to other members of the 
Of?p category, this star does not present large 
photometric changes or impressive variations of 
the optical spectra occurring on very long peri- 
ods (e.g. 538d for HD 191612, -55yrs for HD 108). 
Instead, smaller-amplitude variations of the opti- 
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cal spectrum occur ring on shorter timescales (7 d) 
have been detected ( Naze et al. Il2008a[ Naze et al 



2010l ). The obse rved magnetic field o f HD 148937 
first reported bvlHubrig et al.l (120081) and recently 
confirmed by IWade et al. I (|201l[ ). also seems to 
present no or only very small variations over the 
7d-period. Such a lack of large variations, in con- 
trast with other objects of the same category, does 
not necessarily rule out the Magnetic Oblique Ro- 
tator (MOR) model that is proposed to explain 
the Of?p peculiarities. Small variations could in- 
deed be attr ibuted to a pole-on geometry, as see n 
from Earth (|Naze et al.ll2010l : IWade et al.1l201l[ ). 
though such a geometry would be at odds with the 
appearance of the surrounding ejected nebula. 

In the X-ray domain, the high-energy emission 
of HD 148937 was first detected with Einstein, and 
then re-observed by ROSAT and XMM-Newton 
(jNaze et al. I l2008aL and references therein) . As 



for the two other famous members of the Of?p 
category, HD 108 and HD 191612, the spectrum 
of HD 148937 appeared at the same time over- 
luminous ( log ( L y / ^rqt.) ~ — 6 rather than -7 , 
Naze et al. I l2008at iNazel 120091 : iNaze et ablboilh 
and rather soft (dominant plasma temperature 
at 0.2 or 0.6 keV). These properties are at odds 
with what is expected from both mag netic ob- 
jects such as ^^OriC ( Naze et al.ll2010l and ref- 
erences therein) and typical 0-stars. To go fur- 
ther into the analysis and fully characterize the 
physical properties of HD 148937, high-resolution 
data of high quality are needed. However, the 
high-resolution data from the RGS instrument 
aboard XMM-Newton were very noisy, though 
they seemed to indicate the presence of broad 
X-ray lines (FWHM^2000kms~^), in agreement 
with prediction fr om the 'standard' w ind-shock 
model of 0-stars ( Naze et al. I |2008a). This is 



indeed puzzling since HD 148937 and the other 
Of ?p stars have mag netic confinement p arameters 
as large as 6*^ Ori C (jNaze et al.l l2008b[ ) and one 
would thus expect rather narrow X-ray lines. New 
data were thus essential. 

This paper reports on the analysis of Chandra 
high-resolution data and for completeness we re- 
visit the XMM-Newton RGS data. The observa- 
tions are presented in Sect. 2, the properties of 
the X-ray lines and overall high-energy emission 
are analyzed in Sects. 3 and 4, respectively, while 
Sect. 5 summarizes the results and concludes. 



2. Observations and data reduction 
2.1. Chandra 

HD 148937 was observed by Chandra on June 
5 2010 using the HETG-I-ACIS-S (ObsID=10982, 
PI Y. Naze). The available pipeline data (level 2) 
were not re-processed, as they had been reduced 
with a recent reduction software. They were fur- 
ther reduced with CIAO 4.2 and CALDB 4.3.1. 
Following the recommendations of the Chandra 
teani^, wc first extracted individual first-order 
MEG and HEG spectra for the source and back- 
ground using dmtype2split and tg_bkg, then calcu- 
lated response matrices for both first-orders using 
mkgrmf and fullgarf. The spectra and arf files from 
orders ~1 and +1 were added for each instrument 
using add_grating_orders, and the rmf file for the 
-|-1 order was used as general response matrix for 
each instrument (since the differences between the 
two rmf are very small). For the 0th order data, 
we further defined a source region of 15 pixels ra- 
dius and a surrounding annular background region 
of radii 15 and 45 pixels. We then used the task 
specextract to get the 0th order spectra and its 
corresponding response matrices. 

The total effective exposure was 99ks, and the 
observed count rates in the 0.5-2.5 keV energy 
band are 1.97 ±0.17 x lO^^^tss-i, 5.31 ±0.15 x 
10-2ctss-\ and 4.41 ± 0.07 x lO-^ctss"! for 
HEG, MEG, and 0th order, respectively. In the 
0.5-10. keV range, the count rates (resp. net num- 
ber of counts) are 2.43 ± 0.18 x 10~^ ctss^^ (resp. 
~2400), 5.88 ±0.16 x lO-^ctss"! (resp. -5800), 
and 5.90 ±0.08 x lO^^ctss^i (resp. -5800) for 
HEG, MEG, and 0th order, respectively. It may be 
noted that count rates in the 0.5-10. keV range are 
significantly higher than those in the 0.5-2.5 keV 
energy band, indicating the presence of some hard 
energy flux. This could be due either to a high ab- 
sorption or the presence of very hot plasma, which 
is actually the case here (see below). 

For the line analysis, HEG and MEG spectra 
were re-binned by a factor of 4 and 2, respec- 
tively, to improve the signal-to-noise without de- 
grading the resolution. For the global fitting, the 
HEG, MEG and 0th order spectra were grouped 
to get a minimum o f 20 c ounts per bin. FoUow- 
ing lzhekov fc Fallal ()2007f ). the MEG and HEG 
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spectra near the intercombination and forbidden 
lines in the hehum-hke triplets were re-binned so 
that these two lines fall into one large bin. This 
avoids fitting problems as Xspec thermal mod- 
els (e.g. apec) do not take into account the po- 
tential f /i variations if a strong UV emission is 
present. For both binnings, all available data 
(HAG-hMEG for hne fitting and HEG+MEG-^Oth 
order for global fitting) were fitted simultaneously, 
using a statistics and background subtraction. 

2.2. XMM-Newton 

XMM-Newton observed HD f 48937 on Feb. 
25 200f. Three archival datasets with RGS 
data exist, with ObsID 002 2f40f0f, 002 2f4050f, 
0022f4060f. Compared to iNaze et al. I (floOSa), 
the XMM reduction software has significantly 
improved, especially for the grating data. The 
full observation was therefore downloaded and 
reduced once again with SAS vlO.0.0, following 
the recommendation of the XMM tean|l. The 
data from each instrument and each order were 
combined using the task rgscombine. A high 
background affects all three exposures, but with- 
out a clear flare: the whole exposure (38ks for 
RGSl-orderl and 36ks for RGS2-orderl) was thus 
kept. The overall count rates in the 0.5-2.5 keV 
energy band are 3.71 ± 0.25 x lO^^ctss^^ and 
6.69 ± 0.25 X 10~^ctss~^ for the first orders of 
RGSl and RGS2, respectively. This corresponds 
to net number of counts of about 1400 and 2400 
for the first orders of RGSl and RGS2, respec- 
tively. A grouping of the channels by a factor of 3 
was used, to improve the signal-to-noise without 
degrading the resolution. Note that (1) only first- 
order data were considered since the second-order 
data are very noisy; (2) only RGS2 data are avail- 
able at the wavelength of the neon lines because 
of a CCD failure in RGSl. 

3. Line properties 

Both RGS and MEG spectra are shown in 
Fig. [TJ The strongest lines are NexA12A 
and OVII1AI9A for the RGS and SxvA5A, 
Si XIV A6.2A, Si XIII A6.7A, Mg XII A8.4A, Mg xi A9.2A, 
and Ne X AI2A for the MEG. It may be noted that 
lines of highly-ionized iron are detected, which is 



not common for 'normal' O stars. 

In Fig. [21 we also compare the spectrum 
of HD 148937 to those of the magnetic object 
e^OriC (ObsID 3,4), and two 'normal' O stars, 
15Mon (ObsID 5401,6247) and H P 206267 (Ob- 
sID 18 88,1889), from the sequence of lWalborn et al. 
(|2009l ). Of these four stars, 15Mon suffers 
from the smallest interstellar absorption (E(B- 
V)=0.03), while ^^OriC is more extincted with 
E(B-V)=0.29 and the last two objects display 
the highest extinctions, with E(B-V)^0.5. This 
implies that differences between HD 206267 and 
HD 148937 cannot be due to different interstel- 
lar absorptions, and that the hard character of 
0^ Ori C is not due to a high extinction. 

From Fig. [2l several features are obvious. The 
Ne X line largely dominates over the Ne ix line in 
the spectra of 0^ Ori C and HD 148937. Though 
this pair of lines is the most separated one in wave- 
length, hence potentially the most affected by ab- 
sorption effects, this difference is not due to inter- 
stellar absorption, as HD 148937 and HD 206267 
share similar interstellar extinctions. A similar 
feature is seen for the H-to-He line ratio of mag- 
nesium. Regarding silicon, the Si xiv line is not 
clearly detected in the "normal" 0-stars, contrary 
to Si XIII, while it is of similar strength to Sixiii 
in HD 148937 and even dominates over Si xiii for 
0^ OriC. Finally, the sulfur lines do not appear in 
the spectra of "normal" 0-stars (at the sensitiv- 
ity limit of the data), but are obvious in 0^ OriC 
and HD 148937. It should be noted however that 
S XIV dominates over S xv in HD 148937, whereas 
the reverse situation is seen in 0^ OriC. 

The reversal of the S, Si, Mg and Ne H-to- 
He line ratios between normal stars and mag- 
netic objects is thus clearly seen, the latter hav- 
ing higher ionization. However, it is less ex- 
treme for HD 148937 than in the case of 0^ OriC. 
HD 148937 may nevert heless participa t e in th e di- 
chotomy discussed by IWalborn et al. (I2OO9I) and 
Zhekov fc Paiial (|2007l ). 



3.1. Shifts and widths 

The rest wavelengths of these lines were taken 
from ATOMD^, and Xspec vl2.6.0.q was used 
for the fitting. For He-like triplets (S xv, Si xiii. 
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Mgxi), the shifts and widths of the hnes were 
forced to be equal for the 3 components; whereas 
for the H-hke doublets (Sixiv, Mgxii, Nex, 
Oviii), the relative line intensities were in ad- 
dition tied to their atomic data value at the peak 
emissivity. 

The signal-to-noise of the data prevents us from 
doing detailed line profile modelling: the lines were 
thus fitted by simple gaussians. The fitting results 
are provided in Tableland shown on Figures [3] to 
m For each line, the rest wavelength/energy, shift, 
width, and intensities are tabulated, together with 
their associated Icr error. Note that these errors 
are often asymmetrical: the value shown here al- 
ways is the largest value. The last column provides 
the intensities corrected for an interstellar ab- 
sorption of 4x10^^ cm~^ ( biplas k Savage! Il994l: 



on Chandra data are much smaller than those re- 



Morrison fc McCammonlll983l) . Note that a fit by 
a triangular profile was a l so per formed, for con- 



sistency with lNaze et al. I (|2008al ): the results are 



indistinguishable from those reported in Table [TJ 
The results for the Chandra data are rather 
stable from line to line: the average shift is 
49 ± 32kms"\ and the average FWHM is 
827±90kms"i (note that quoted errors are er- 
rors on the mean, not dispersions around the 
mean). The lines are thus not significantly 
blue- and only slightly red-shifted , in agree- 
ment w ith the null shift found by Naze et al. I 
(|2008al ) . Thi s is not surprising since, apart from a 
few objects (iKahn et al.ll200lt iRauw et al.] 12002 : 



Cohen et al.l 120111 ). 'norm al' Q-stars generally 



show unshifted lines (see Waldron fc Cassinellil 



20071: iGiidel fc Na^l2009l and references therein). 



The observed widths can also be compared to val- 
ues reported for mas sive stars in the literature (see 
Giidel fc Na^ l2009l and references therein): the 



first 0-stars observed at high-resolution yielded 
mixed results, some having broad lines (FWHM of 
2000kms-i for CPup and ITOOkms"^ for COri), 
other narrower lines (FWHM of 860km s"-'^ for 
JOri, 800kms~"'^ for COph); while global studies 
later showed that FWHM< Voo for most stars. 
Indeed, the average width found here is below the 
wind terminal velocity (see Table [2]) , in line with 
the "narrow" cases mentioned above, but it is not 
as extreme as the FWHM=600kms~^ repo rted 



erratum m 



Gagne et al 



for most li nes of 6^ Ori C ( Gagne et al.ll2005l and 
l2005l) 



ported for the same star by iNaze et al. I ()2008a , 
about 2000km s-^) using XMM-Newton RGS 
data. Since the reduction software has signif- 
icantly changed in the last years, we have re- 
reduced the RGS data and re-fit the lines (see 
Sect 2.2, Figs. [5]fc|6]and Table [J): once again, 
the oxygen line appears resolved and we found 
large values for the widths. The differences thus 
appear real, but two remarks must be made: (1) 
the signal-to-noise of the XMM data is much lower 
than that of the Chandra data, which could blur 
the results (indeed, the error bars are much larger, 
see Fig. |6l), and (2) the O Vili line is generally ex- 
cited by lower-temperature plasma. As the X-ray 
emission is most probably not isothermal (see next 
section), the data difference could simply reflect 
the presence of s everal plasrnas. Indeed, in the 
case of 0^ OiiC ( Gagne et al. 2005t Gagne et al 



20051 ). the Oviii width was found to be larger 
than that of shorter wavelength lines by a factor 
of 2.5 (which is 40% higher than what we found for 
HD 148937). Therefore, what the Chandra data 
tell us is that the high-temperature plasma, most 
probably linked to magnetic confinement close to 
the star, displays narrow lines, i n agreement wit h 
theoretical predictions (see e.g. ud-Doula 2008[) . 
This solves at least one problem in the under- 
standing of the X-ray emission from Of?p stars. 

3.2. Ratios 

Using the line intensities, we calculated the 
= f Ih Q = if + T^jf, and H-to-He-like ra- 
tios, which we compare to theoretical values for 
deriving plasma location and temperature. 

We extracted the theoretical emissivity in- 
formation for the H-like and He-like lines from 
ATOMDB v2.0.0 (iFoster et al. Il201l[ ). Note that 



However, the line widths measured for HD 148937 



ATOMDB calculations are done considering col- 
lisional ionization equilibrium and that the used 
emissivities are not limited to those of the lines 
themselves as several other lines blend with the 
doublets and triplets. For our calculations, we 
have therefore summed the emissivities of all 
lines found over an interval Aq ± 0.5 x resolution 
around the lines of interest, using resolutions of 
1.2x10-^1, 2.4xlO-2A, and GxlO'^A for HEG, 
MEG and RGSl-orderl, respectively, which corre- 
spond to the FWHMs of the instrumental broad- 
ening. 
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The observed H-to-He-like ratios of Si and 
Mg favor teniperatures log(r) of 7.08±0.02 and 
6.92±0.05, respectively. These high temperatures 
(12 and 8MK) are l arger than those of "n o rmal" 
OB stars quoted in IWaldron fc CassineU? (200 
see also the erratum in IWaldron & Cassincll: 



1 

I 



20081 ). but similar to those reported for the col- 



liding wind binary Cyg 0B2 #8a and the mag- 
netic star rSco in the same reference. It should 
be noted, however, that 9^ Ori C displays much 
higher tempe ratures (10-12MK for Mg an d 15MK 
for Si). As iWaldron fc Cassinellil (|2007t ) noted, 
those temperatures agree well with the (ranges of) 
temperatures of maximum emissivities of these H- 
and He- like ions (this is also true for O, see below). 

The g value that we derive for HD 148937 in- 
dicates similar log(r) for the different elements: 
6.95±0.25, 6.82±0.2 and 6.87±0.5 for S, Si and 
Mg, respectively. While the G ratios for Mg and 
Si are similar to thos e of OB stars (w hatever their 
nature) reported by IWaldron fc Cassi ncUi (200^ 



20081 ) ■ the G ratio derived for S can only be com- 



pared to those of the magnetic objects t Sco and 
0^OriC, as well as to those of giant and main- 
sequence stars (only supergiants having higher ra- 
tios). It may be noted that the above tempera- 
tures corresponds to ~0.7keV, which agrees well 
with the "cool" temperature of the "hot" global 
model (see next section). 

At temperatures log(r) of 6.8-7.0, TZq (cal- 
culated using ATOMDB as above) is about 2 
for S, 2.65 for Si, and 2.42 for Mg. It is well 



known that TZ 



^0 —4 



7+riUr, 



where 



we can neglect the density dependence for mas- 
sive stars. The UV flux depends on the stel- 
lar output but also the dilution factor. We de- 
rive the UV flux for HD 148937 using the model 
(Te//=40kK and log(g)=4.0) in the grid of O- 
star atmosphere models calculated with CMFGEN 
(iHillier fc Milleil lT998'R which is closest to the pa- 
rameters derived from a d edicated atm osphere fit- 
ting of HD 148937 (Naze et al. ll2008aL and Table 
[2|). After averaging the flux in the velocity in- 
terval [-2000,0]kms~-'^near the rest wavelengths of 
the 2^ Si — 2'^ Pi, 2 transitions ( follow ing a sim- 
ilar path as iLeutenegger et aL 1 liooi), see also 
Table|4]), we then derived a formation radius of 
1.9±0.4i?* from the 7^ ratios of Si, and <3.2i?* 



for Mg (the ratio itself corresponding to 1.7i?*j3. 
The sulphur lines are too noisy to provide a mean- 
ingful constraint on the formation radius. Such 
rather close radii are similar to thos e generally 
found for 0-stars (ICiidel fc Na ze"2009^. including 
6'^OriC (iGaene et al l 120051: iGagnc et al. 200^: 
therefore, it cannot be used to discriminate be- 
tween various X-ray emission models. It may be 
worth noting, though, that our values fit well the 
temperature vs forma t ion radius r elation found by 
Waldron fc Cassinellil (|2007ll2008l ). 



Following Hucne moerder et al. ( 20091) , we de- 
rive an abundance ratio Mg/Si, relat i ve to the 
solar abundance of Anders fc G reves"s3 ( 19891 ) , of 
0.74±0.11 using the H-like and He-like resonance 
lines. Non-solar abundances in CNO elements are 
not entirely surprising since HD 148937 displays 
both a nitrogen ove rabundance and an enriched 



surrounding nebula (jNaze et al. I l2008al . and ref- 
erences therein). However, changes in the Mg/Si 
ratio are not expected (see also next section). 

Finally, we may try to use the ratio of H-to- 
He lines for oxygen to check whether our hy- 
pothesis of a cooler plasma is sensible. The 
Ovii triplet is very noisy in the RGS data, 
but the fiux in the r line can be constrained to 
2.4±0.6xl0-5phcm-2s-i (observed) or 69±19x 
10~^ phcm^^ s~^ (corrected for absorption) . This 
corresponds to a H-to-He ratio of 0.67±0.20, yield- 
ing a temperature of log(r) ^ 6.35. This is much 
lower than the temperatures found above: it would 
be difficult to reproduce this measurement in the 
context of an isothermal plasma, thereby confirm- 
ing our hypothesis of a multi-temperature plasma. 
This result is certainly far from perfect because of 
its high uncertainty, but it constitutes a hint to- 
wards the non-uniform temperature of the X-ray 
emitting regions. 

4. Global fits 

Global fits were simultaneously made on the 
HEG, MEG and 0th order spectra, using a bin- 
ning ensuring at least 20 counts per bin (see Sect. 
2). Results of the fits are presented in Table [5] 



* http://kookaburra.phyast.pitt.odu/hillier/web/CMFGEN. htm 



^Quoted errors take into account only the statistical errors 
of the measured fluxes of the f and i lines (and considering 
error propagation), not errors on the predicted CMFGEN 
EUV fluxes, which are difficult to quantify without actual 
observations of HD 148937 at these EUV wavelengths. 
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Note again that the quoted la errors are some- 
times asymetricah the value shown here always is 
the largest value. 

We used two sets of models that both assume 
thermal plasma in coUisional ionization equlibrium 
(CIE). First, we used models with discrete temper- 
ature components {apec+apec or apec+apec+apec 
in Xspec). The second set models the emis- 
sion of a distribution of radiative shocks as in 



Zhekov fc Pallal (|200iT ). 

These emission models were multiplied by 
two absorption components: one corresponds to 
the fixed ISM absorption (model w abs, — 



Diplas fc Savagd 119941) ■ the other 



to potential additional local absorption (model 
phabs). Considering the results from the line 
fitting, the resulting absorption-|-emission spec- 
trum was convolved with a constant-velocity 
Gaussian broadening using the gsmooth model 
in Xspec (width at 6keV=7.024eV and a{E) = 
a{&keV) X E/Q) while keeping the line shift to 
zero. 

In each case, we first perform a fitting using 
solar abundances (with a dopted abundances from 
Anders fc Grevesse Il989^ and we then varied the 
Ne, Mg, Si, S, and Fe abundances (i.e. those from 
elements displaying strong lines) to improve the 
fitting quality. Ne, Mg, and Si remain close to 
solar values, while S appears overabundant com- 
pared to solar and Fe subsolar. The major fit 
improvement in this case mostly comes from the 
changing iron abundance - the improvement from 
other elements is not formally significant. 



As found bv lNaze et al. I (|2008af ). the presence 
of hot plasma is needed for a good representation 
of the X-ray emission from HD 148937. However, 
we found again the same "duality" in the 2T fits, 
with both a "cool" model (0.3keV-|-2keV) and a 
"hot" model (0.6keV-|-2.3keV) giving rather sim- 
ilar results. A three-temperature model rather 
yields best-fit plasma at 0.25, 0.8, and 2.4 keV, 
but the improvement in , compared to 2T mod- 
els, is not very large. The relative importance 
of the two temperatures differs in these models, 
though: the 0.3 keV component is largely domi- 
nant in the "cool" and 3T models, whereas the 
two temperatures of the "hot" model and the two 
higher temperatures of the 3T model have simi- 
lar emission measures. As usual, there is clearly 
an interplay between temperatures, normalization 



factors and absor ptions, but compared to "nor- 
mal" 0-stars (e.g. lNazeir2009h . the strength of the 
hottest component is much higher here. 

The distributed shock model favor similar peak 
temperatures (0.25, 1. and 3keV , Fig. [7t, again 
in contrast to "normal" 0-stars (Zhekov fc Palla 



20071 ). but with more contributions from the high- 



est temperatures. This apparent difference in the 
shape of the distribution of emission measure of 
CIE plasma and that of radiative shocks deserves 
some comments. As described in IZhekov fc Palla 
((20071) , at each temperature value the model with 
distribution of radiative shocks presents the to- 
tal amount of emission measure of shocks with 
that given postshock temperature. Since the post- 
shock region in a radiative shock is a temperature- 
stratified region by itself (occupied by plasma with 
temperature lower than the immediate postshock 
temperature, which is the one represented on Fig. 
[71 explaining the skewing towards higher temper- 
atures), then the actual distribution of thermal 
plasma in the entire X-ray emitting region is dif- 
ferent from that of the radiative shocks. For such 



a reason, iZhekov fc Palla (200Z) emphasized that 
it is not correct to use canonical models of dis- 
tribution of emission measure of CIE plasma to 
deduce the possible immediate postshock temper- 
atures (and shock velocities) in the analysis of the 
X-ray emission from hot massive stars: the shocks 
in their m assive winds are expec ted to be radiative 
(see § 2 in lZhekov fc Pallall2007 and the references 
therein). 

Using the results from the global spectra fits, 
we can calculate the so-called total stellar 'cloudi- 
ne ss' for HP 148937 (se e § 5 and Appendix A 
in IZhekov fc Pallal 120071 ) . Interestingly, we get 



a value of 0.008 which together with the zero- 
velocity line shift for its X-ray spectral lines (see 
previous section) puts HD 148937 exactly next to 
0^ Ori C on the corresponding plo t for a sample of 
hot m assive stars (see Fig. 4 in IZhekov fc Palla 
20071 ). Indeed, when comparing the radiative 



shock model of HP 1 48937 with the sample of 
iZhekov fc Pailal (|2007l ). it is obvious that the only 
object comparable to HD 148937 is 0^ Ori C (and 
maybe rSco). However, whilst the temperature 
distributions of both objects for a radiative shock 
model show peaks at high temperatures {kT > 
1 keV), HD 148937 displays in addition a consider- 
able amount of shock emission at temperatures be- 
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low 1 keV - the emission below 1 keV is rather neg- 
ligible for 9^ Ori C (c ompare Fig. [7] and Fig. 2 of 
Zhekov fc Pallall2007[ ). Similarly, the temperature 
distribution of CIE models show a smaller impor- 
tance of the high-temperature peak compared to 
the low-temperature one for HP 148937 (compar e 
Table [5] here and Table 2 in iGagne etaP l2005[ ) . 
Indeed, a 2T fit of 6*1 Ori C yields values of 0.7 
and 2.8 keV with the hottest component having 
an emission measure 4 times larger than that of 
the coolest one, whereas the emission measures are 
similar for the "hot" model of HD 148937, which 
has temperatures similar to those of the 2T fit of 
0^ Ori C. This suggests that most of the overlu- 
minosity of HD 148937 comes from relatively soft 
plasma. 

In general, the absorption does not seem to 
exceed the ISM value by much. The derived 
abundance values are in general consistent with 
solar abundances, except for iron (but the so- 
lar abundance of i ron has been revised since 
1989t) : the global fits con- 



Anders fc Grevesse 



firm the 3/4 solar abundance of Mg/Si found 
from the line ratios. Finally, the flux observed by 
Chandra is consistent with that derived from the 



XMM -Newton low-resolution spectra (|Naze et al. I 
2008al ). confirming the lack of large variations in 
the X-ray domain. 

5. Summary and conclusion 

HD 148937 is a peculiar object belonging to the 
intriguing category of Of?p stars. The presence of 
a magnetic field in these stars, together with their 
optical variations are generally explained within 
the framework of the Magnetic Oblique Rotator 
model. In 2008, a first X-ray investigation of 
HD 148937 using XMM-Newton had revealed pe- 
culiarities in its X-ray spectrum too. Unfortu- 
nately, no detailed, high-resolution data of good 
quality was available until we obtained new data 
with Chandra-HETGS. 

These new data enable us to study HD 148937 
in more detail. The X-ray lines detected by Chan- 
dra display an average shift of 49 ± 32kms~^, 
and an average FWHM of 827 ± 90kms-i. The 
lines from high-Z elements, whose emissivity peaks 
at high temperatures, thus appear rather narrow, 
suggesting that they could be linked to magnet- 
ically confined winds. If the larger width of the 



O VIII line found from the noisier RGS data is real, 
it means that there is no unique source or location 
for the hot plasma. This is a similar situa tion to 
what was found for 6*1 OriC (ICaene et all [20051) - 
though the FWHM of the high-Z lines was even 
narrower for this star, about 600 km s""'^. 

While RGS data already provided hints of a 
suppressed forbidd en line for the Neix triplet 
(|Naze et al. Il2008a[) . the good quality HETG spec- 
tra confirm these hints and provide solid evidence 
for suppressed forbidden lines in two other He-like 
triplets, those of Sixiii and Mgxii. These He-like 
triplets of Mg and Si indicate a formation radius 
of about one stellar radius from the photosphere, 
which is rather close to the star's surface, but 
not uncommon in "n ormal" or magnetic 0-stars 
(jGiidel fc Nazell2009^ . 

The temperatures derived from H-to-He like 
ratios are higher than for "normal" 0-stars, 
but still not as extreme as th ose of 6^ Ori C 
(|Waldron fc Cassinelli|[2007ll2008l) . This indicates 
the presence of hot plasma, which is confirmed in 
global fits. Global fits further show that this very 
hot component is much stronger than for "nor- 
mal" 0-stars but still much less dominant than 



in the extreme case of Ori C ([Zhekov fc Palla 
2007t iGagne et allbOOSl) . The comparison of H- 



like and He-like line pairs supports this view, as 
HD 148937 displays a higher ionization than "nor- 
mal" O stars, albeit not qui te as extreme as in 
6*1 OriC (IWalborn et al.ll2009l) . 

The fitting also shows that there is no large 
additional absorption over the ISM value, and 
that the overall flux agrees within 10% with that 
recorded 8.5 years before by XMM-Newton. Re- 
garding abundances of Ne, Mg, Si, and S, the 
global fits favor close-to-solar abundances (with a 
slight excess in sulphur and a deficit in iron). 

In conclusion, the new, good quality high- 
resolution data obtained by Chandra yielded addi- 
tional clues in favor of the MOR model: relatively 
narrow X-ray lines (FWHM~800kms~^) and a 
si milar 'cloudiness ' for th e radiative shock model 
of IZhekov fc Fallal (|2007h . plus a confirmation of 
the presence of very hot plasma. However, abun- 
dant soft plasma is also present, in contrast with 
the MOR-prototype 9^ Ori C: in every respect (e.g. 
line widths of high-Z elements or temperature dis- 
tributions), HD 148937 thus appears less extreme 
than 9^ OriC, despite their similar magnetic con- 
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finement parameters. Further MHD modelling of 
confined winds are definitely needed to explain 
this difference. 
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Fig. 1.— The MEG and RGS spectra with the 
strongest lines identified. Note that (1) the wave- 
length range is not identical in both panels be- 
cause of the different sensitivities of the consid- 
ered instruments and (2) the spectra have been 
background subtracted. 




X (A) 



Fig. 2. — MEG spectra of a selection of OB stars 
(spectral channels were binned by a fac tor of ten 
for all stars, as in IWalborn et all l2009l ). Actual 
fluxes have been normalized so that the strongest 
lines peak to 1, and the most important H- and 
He-like pairs are labelled. 
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Fig. 3. — The observed H-like line profiles and 
their best fit models. Dotted lines are drawed at 
velocities of ~Voo, and -\-Voo- 
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Fig. 7. — The distribution of emission measure as a function of immediate postschock temperature for a 
model with distributed radiative shocks. 
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Fig. 4. — The observed He-like line triplets and 
their best fit models. Dotted vertical lines indicate 
the rest wavelengths of the rif triplets, and the 
small horizontal line below the instrument's name 
has a length of 1000km s~^. 



Fig. 5. — The line profiles observed with RGS and 
their best fit models. The left panel shows the 
RGS2 data for Ne X while the right panel displays 
both RGS for O viii, the thin black and thick red 
line corresponding to RGS 1 and 2, respectively. 
Dotted lines are drawed at velocities of — Woo, 
and +Voo- 
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Fig. 6.— The best-fit FWHMs and line shifts. 
Fihed circles refer to Chandra data, open circles to 
RGS; horizontal lines indicate the average values 
and their associated ±lcr interval. 
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Table 1: Best- fit properties of the observed lines. 



Ion 


A 


E 


shift 


FWHM 


Obs. Intensity Cor. Intensity 




(A) 


(kcV) 


(km s^"'") 


(kms-i) 


(phcm^^s^^) (phcm^^s^^) 














Sxv 


5.03873 


2.46062 


88±225 


962±809 


(9.3±2.8)xl0"'' (10.2±3.1)xl0"'' 




5.06649 


2.44714 






(2.6±2.6)xl0"'^ (2.9±2.9)xl0~® 




5.10150 


2.43035 






(5.5±2.6)xl0-'^ (6.1±2.9)xl0-6 


Si XIV 


6.18044 


2.00607 


0±64 


805±252 


(9.30±1.08)xl0"^ (11.0±1.28)xl0"^ 




6.18585 


2.00432 






previous x . 509 


Si XIII 


6.64795 


1.86500 


117±59 


859±160 


(1.46±0.18)xl0-5 (1.80±0.22)xl0"''^ 




6.68819 


1.85378 






(7.52±1.43) x 10^^ (9.30±1.77) x 10"® 




6.74029 


1.83945 






(6.02±1.30) x 10-6 (7.32±1.58) x IQ-^ 


Mgxii 


8.41921 


1.47263 


33±57 


769±152 


(8.78±1.15)xl0-6 (12.6±1.65)xl0~6 




8.42462 


1.47169 






previous X 0.507 


Mgxi 


9.16875 


1.35225 


42±123 


894±297 


(1.01±0.25)xl0-5 (1.59±0.39)xl0-5 




9.23121 


1.34310 






(7.0±2.6)xl0-6 (ll.l±4.1)xlO-6 




9.31434 


1.33111 






(0.5±1.4)xl0-'^ (0.8±2.3)xl0-6 


Nex 


12.13210 


1.02195 


-15±150 


876±275 


(1.73±0.33) x 10-5 (4.31±0.82) x lO^^ 




12.13750 


1.02150 






previous X 0.502 


XMM 












Nex 


12.13210 


1.02195 


-18±510 


1815±1124 


(4.0±0.7)xl0-5 (10.0±1.7)xl0-5 




12.13750 


1.02150 






previous X 0.502 


viii 


18.96710 


0.65368 


-440±180 


1537±505 


(2.9±0.4)xl0-5 (30.7±4.2)xl0-5 




18.97250 


0.65349 






previous X 0.501 



Table 2: Stellar parameters (from iNaze et al. 



2008al ). Note that the wind parameters remain 



uncertain. 



Parameter 


Value 




41±2kK 


logg 


4.0±0.1 


Voo 


-2600km s-i 


log(Af/V7) 


-6 



Table 3: Observed line ratios. 



Element 


R 


G 


H-to-He 


H-to-He (r only) 


S 


2.1±2.3 


0.88±0.48 






Si 


0.79±0.23 


0.92±0.17 


0.48±0.07 


0.92±0.16 


Mg 


0.07±0.21 


0.75±0.35 


0.68±0.17 


1.19±0.33 
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Table 4: Rest wavelengths of the 2 -J> 2 ^p^ ^ 
transitions in ATOMDB v 2.0.0. 



Ion 


Ai(A, 2=^51 ^ 2^Pi) 


A2(A, 2^ Si 2^P2) 


Sxv 


738.3 


673.4 


Si XIII 


865.1 


814.7 


Mgxi 


1034.3 


997.5 



Table 5: Results of the fitting. Abundances are 
in numbers, relative to hydrogen and relative to 
solar. 



Parameter 






2T 


cool 


2Th 


ot 


3T 




Rad. shocks 








(solar) 


(var.) 


(solar) 


(var.) 


(solar) 




(var.) 




Nji (10^^ cm-^) 




0.34±0.03 


0.14±0.04 


0.13±0.02 


O.iO.Ol 


0.27±0.03 


0, 


.13±0.06 


0.11±0.07 


kTi (kcV) 






0.28±0.01 


0.37±0.02 


0.64±0.01 


0.65±0.01 


0.24±0.02 


0, 


.30±0.03 




normi (10~ 


^ cm" 




10.9±1.47 


4.47±0.90 


1.89±0.16 


1.78±0.16 


7.28±2.30 


2, 


,89±1.25 


8.25 


m (keV) 






1.91±0.04 


2.04±0.06 


2.34±0.06 


2.44±0.08 


0.78±0.03 


0, 


.79±0.04 




norm,2 (10~ 


^ cm" 




2.66±0.06 


2.39±0.09 


1.93±0.05 


1.83±0.07 


1.53±0.15 


1, 


.31±0.19 




m (keV) 














2.33±0.08 


2, 


,52±0.10 




nor 1713 (10~ 


^ cm" 












1.91±0.07 


1, 


.72±0.08 




Ne 








0.92±0.12 




0.96±0.14 




0, 


.90±0.15 


0.86±0.16 


Mg 








0.84±0.15 




0.86±0.10 




0, 


.99±0.12 


0.85±0.12 


Si 








1.15±0.11 




1.07±0.11 




1, 


.15±0.12 


1.03±0.11 


S 








1.41±G.19 




1.73±0.22 




1, 


.82±0.24 


1.86±0.25 


Fe 








0.53±0.06 




0.59±0.06 




0, 


.75±0.09 


0.71±0.11 


Flux (10-12 


'erg cm ^s i) 


2.70 


2.70 


2.73 


2.76 


2.77 




2.79 


2.99 


(dof) 






1.02 (556) 


0.94 (551) 


0.84 (556) 


0.74 (551) 


0.73 (554) 


0. 


,69 (549) 


0.61 (534) 
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